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Abstract 

We show that all correct results obtained by applying quantum field theory to neutrino 
oscillations can be understood in terms of a single oscillation formula. In particular, the 
model proposed by Grimus and Stockinger is shown to be a subcase of the model proposed 
by Giunti, Kim and Lee, while the new oscillation formulas proposed by Ioannisian and 
qq ' Pilaftsis and by Shtanov are disproved. We derive an oscillation formula without making 

\& • any relativistic assumption and taking into account the dispersion, so that the result is valid 

for both neutrinos and mesons. This unification gives a stronger phcnomcnological basis to 
the neutrino oscillation formula. We also prove that the coherence length can be increased 
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without bound by more accurate energy measurements. Finally, we insist on the wave packet 
interpretation of the quantum field treatment of oscillations. 



Q*. 1 Introduction 

Nowadays the study of neutrino oscillations is the particle physicists' best hope to learn how to 
extend the Standard Model, since it provides information on the neutrino masses and mixings 
which are thought to originate in new physics. Neutrino oscillations are an example of flavor 



oscillations, which arise when the particles produced and detected in an experiment are superpo- 
sitions of different mass eigenstates. The interference between the propagating mass eigenstates 
leads to a spatial oscillation of the detection probability of a neutrino with a given flavor. The 
standard oscillation formula, expressing the probability to detect a neutrino of momentum po 
emitted with a flavor a as having a flavor (3 after a propagation on a distance L, reads Q 

-i \ pQ 3 Lj , (1) 

where rriij are the neutrino mass eigenstates and V is the matrix relating the flavor to the mass 
fields. This phenomenon can explain the depletion of the solar electron-neutrino flux ||, as well 
as the up/down asymmetry of the neutrino flux produced by cosmic rays in the atmosphere Q. 
Moreover, the long baseline experiment K2K hints to a muon-neutrino depletion on a terrestrial 
distance Q , while controversial evidence for neutrino oscillations has appeared in the laboratory 
experiment LSND H. 
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In view of the importance of neutrino oscillations, it is rather unsettling that the neutrino os- 
cillation formula ([!]) used to analyze experimental results is nearly always derived in a completely 
inconsistent way. In particular, this derivation usually starts with unlocalized plane wave states 
and ends, strangely, with a space-dependent oscillation formula. This hocus-pocus requires at 
least two arbitrary assumptions: the equal time prescription, stating that the propagation time 
is the same for the different mass eigenstates, and the classical propagation condition, x = vt 
(see H for a review). Additional assumptions have generated endless discussions about the 
equality of the energy-momenta of the different mass eigenstates. Moreover, the rejection of the 
equal time prescription by some authors has lead them to predict oscillation lengths differing 
from the standard result by some factor (a factor 2 is often quoted). 

A more consistent derivation of the neutrino oscillation formula has been done in the in- 
termediate wave packet model, in which the neutrino mass eigenstates are represented by wave 
packets propagating in space-time J?]] . Although the classical propagation condition can be dis- 
pensed with, this method implicitly requires the use of the arbitrary equal time prescription in 
order to obtain the standard result for the oscillation length. Moreover, this derivation is not 
wholly consistent because it involves flavor eigenstates which cannot be well-defined §§. Finally, 
while the 3-momentum uncertainty is naturally included in this model, the energy uncertainty 
has to be put in by hand, with the result that there is no agreement on whether the coherence 
length has an upper bound or not || [l^] . 

In order to solve the various problems of the intermediate wave packet model, different 
models using quantum field theory have been proposed. They can be grouped in four categories. 
The first category groups the external wave packet models, the best example of which is given by 
the seminal paper by Giunti, Kim, Lee and Lee fll|| . In these models, the neutrino is considered 
as an internal line of a Feynman diagram, propagating between a source and a detector, which 



are represented by in- and outgoing wave packets [11, 12, 13, Hj. A variation on this theme is 
the Kiers- Weiss model, in which the external wave packets are replaced by quantum oscillators 
[Oj. The second category groups the stationary boundary conditions models, the best example 



of which is given by the Grimus-Stockinger model [16]. They are very similar to the external 
wave packet models, except that the wave packets are replaced by states with a well-defined 
energy, i.e. stationary states jy], [u], [L7], |l^, |i~9|| . The third category groups models where the 
neutrino is represented by its propagator and coupled at production with a source, though not 



coupled at detection with anything [20, 21]. Note that this kind of model requires the equal 
time prescription, as in the intermediate wave packet model, in order to obtain the standard 
oscillation length, since the detection mechanism is left unspecified. The fourth category includes 
the Blasone-Vitiello model, in which the construction of a Fock space of flavor eigenstates is 
attempted ||. 

While few authors deny that the most rigorous treatment of oscillations is done in a quan- 
tum field framework, the variety of such models has not made a good case for their widespread 
acceptance. Besides, these models sometimes yield conflicting oscillation formulas, either be- 
cause of their different assumptions or because of their different approximation methods. The 
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main aim of this article is to show which of these results are correct and how they can all be 
understood within the external wave packet model. In particular, the Grimus-Stockinger os- 
cillation formula will be seen to be a limiting case of the Giunti-Kim-Lee formula, while the 
Ioannisian-Pilaftsis and Shtanov new oscillation formulas will be shown to be false. At the same 
time, the role of accurate energy measurements in the coherence length will be clarified. 

Contrary to the existing literature, our calculations take into account the dispersion and 
do not resort to a relativistic limit. In this way, the final oscillation formula will be valid as 
well for neutrinos as for K and B mesons (it will be shown that the spin can be neglected). 
This unification of the neutrino and meson oscillation formulas makes it much more difficult 
to accept nonstandard results, since the meson oscillations are well studied and the meson 
oscillation length has been cross-checked with different methods. Another of our concerns will 
be to underline the quantum mechanical interpretation of the external wave packet model. The 
quantum field theory models are often perceived as complicated and obscure, whereas they also 
exist in simple versions and have a clear physical interpretation in direct correspondence 
with intermediate wave packet models. 

The outline of the article is as follows. In the second section, the gaussian external wave 
packet model is reviewed and all useful definitions are given. The flavor-mixing amplitude is 
expressed as a convolution of the neutrino propagator with a function depending on the overlap 
of the external wave packets. This section does not contain new material and follows the notation 



of [12]. In the third section, a new way to evaluate the convolution is explained, leading to an 
explicit oscillation formula. It is then shown that the coherence length, beyond which oscillations 
vanish, can be increased without bound by more accurate energy measurements. As a corollary, 
it is proved that the stationary boundary conditions model is a special case of the external 
wave packet model. In the fourth section, another new method to evaluate the convolution is 
presented, yielding the explicit dependence of the amplitude on time and distance. This method 
allows a clear physical interpretation of the amplitude in terms of wave packets. In particular, 
the propagation range is subdivided into three regimes by two thresholds marking the onset 
of the transversal and longitudinal dispersion of the wave packet. With this method, the new 



oscillation formulas proposed in Refs. ITq] and |21] can be disproved. Finally, we show that the 



Blasone-Vitiello oscillation formula [22] is phenomenologically equivalent to the standard result. 



2 The external wave packet model 

The propagating process of a particle between a source and a detector (indicated by dotted 
circles) is symbolized by the following diagram: 
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Pi(q) — 
Dj(q') 




FIG. 1. Propagation of a particle v between a source and a detector, centered in x P and x D . 

The arrows indicate the momentum flow. Pj represents the set of incoming particles, of 
total momentum q, arriving in the production region (or source), which is centered around the 
point (tp,x P ). Pp represents the set of outgoing particles, of total momentum k, coming from 
the production region, with the exception of the intermediate neutrino v whose propagation is 
studied. Dj, Dp and (t D ,x D ) are defined similarly, but apply to the detection process. The 
interaction points at production and detection are noted x and x', respectively. All particles 
are assumed to be stable. The internal line could also represent an antineutrino, but we shall 
assume that the experimental conditions are such that a quasi-real neutrino propagates on a 
macroscopic distance, transferring positive energy from x to x' . If x'° — x° is a macroscopic time, 
it will be seen that the neutrino v contributes to the propagation, but not its antiparticle. 

Since known neutrinos are much lighter than the experimental energy thresholds, their 
masses can be neglected in the interaction vertices and in the numerator of the neutrino propa- 
gator. In that case, the spin structure factorizes from the sum over the mass eigenstates, so that 
the computation of the oscillation formula can be carried out as if the neutrinos were scalar. If 
one wants to consider nonrelativistic neutrinos, factorization also occurs for nearly degenerate 
masses. The only case where factorization does not occur is nonrelativistic neutrinos with very 
different masses. However, spatial oscillations vanish in these conditions since it is experimen- 
tally possible to determine which mass eigenstate contributes to the amplitude. In that case, 
a treatment with wave packets is not necessary, since the flavor-changing probability does not 
depend on the propagation distance. In the following we shall assume that the neutrinos are 
either relativistic or nearly degenerate in mass. Thus the spin structure can be absorbed in an 
overall multiplying factor. 

The external wave packets are expressed in terms of their momentum-space wave function: 
\tp >= f[dk\i/j(\i, K) |k >, where [dk] is the integration measure in 3-momentum space and K 
is the average 3-momentum. If ^(k, K) is the momentum space wave function of a wave packet 
centered in x = at time t = 0, then 

*(k, K, x , t ) = V<(k, K) e ^(k)to-*k.xo ; 

is the momentum-space wave function of a wave packet centered in xo at time to- 

Without loss of generality, let us choose to work with only one particle in Pi(q), in Pp(k), 
in Dj(q') and in Dp(k'). The extension to a larger number is straightforward and would only 
complicate the notation. The wave packets are built such that those involved in the production 
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of the v are centered at x P at time t P , while those involved in the detection of the v are centered 
at x D at time t D . They are noted 

\Pi> = j[dq}^p I (ci,Q,x P ,t P )\P I (ci)> 

\Pf> = J[dk]*p F (k,K, Xp ,t P )\P F (k)> 

\Di> = J[d< l '}y Dl (< l ',Q > ,x D ,t D )\B I (q > )> 

\D F > = Jldk'}y DF (k',K',x D ,t D )\D F (k , )> . 

Let us first suppose that there is only one mass eigenstate rrij associated to the propagating 
neutrino of FIG. 1. The general formula of the connected amplitude corresponding to this 
process is given by 

A j =<P F ,D F \f (expf -i j d 4 xHi ] J -1\P I ,D I >, 



where Hi is the interaction Lagrangian for the neutrino and T is the time ordering operator. 
Let g be the coupling constant of v with the other fields. Expanding the amplitude to order g 2 , 
and inserting the wave packets expressions, one obtains 



' (27r) 4 p 2 — m A - + te 



where the average propagation time T is defined by T = x® D — x p and the average propagation 
distance by L = x D - x P . 

The overlap function V'liAp) represents the overlap of the incoming and outgoing wave 
packets, both at the source and at the detector. It is defined by 



& ikx 



^(p°,p) = J d 4 x^l d'x'e-^' J[dq\^c l M)e- iqx J[dk]r Po JKK) 

x j [dq']^ in (q',Q')e-^ V J [dk'^J*' e ik '*' M P (q,k) M D (q' ,k') , (3) 
where Mp(q,k) and M£>(q',k') are the interaction amplitudes at production and detection. All 



external particles are on-shell, i.e. q° = Ep ; (q) = w q 2 + m P] , and so on. Note that the overlap 
function is independent of x P and x D . The integrals over x and x' in Eq. @ yield delta functions, 
which impose energy-momentum conservation at the source and the detector. The amplitude 
corresponding to the propagation of an antineutrino from x P to x D can be obtained by the 
substitution p — ► —p. 

Let us now consider the case where the intermediate state in FIG. 1 is a superposition 
of different mass eigenstates. The fields in the flavor and mass bases are related by a unitary 
matrix V Q 
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where Greek and Latin indices refer to the flavor and to the mass basis, respectively. We 
assume that it is possible to tag the initial flavor a and the final flavor j3 of the neutrino. The 
flavor-mixing amplitude A a —,/3(T, L) is defined as the transition amplitude corresponding to the 
propagation over a time T and a distance L of a neutrino initially of flavor a, but detected with 
a flavor (3. 

The flavor-mixing amplitude can be expressed as a linear combination of amplitudes Aj 
corresponding to the propagation of different mass eigenstates: 

A a ^(T,L)=Y / V^A j V ja , (4) 

3 

where Aj is given by Eq. (|2|). 

As the propagation time is not measured in experiments, the transition probability has to 
be averaged over T: 

V a ^(L) ~ Y^ViaVfc V* a V jP j dTAiA* • (5) 

i,3 



A discussion of the proportionality factor can be found in Ref. 13]. 

The next task consists in performing the momentum integration in the amplitude (|2|). This 
cannot be done without doing some assumptions on the overlap function (|3|) , which should be as 
weak as possible. A fairly general oscillation formula can be obtained by approximating the in- 



and outgoing particles with gaussian wave packets |llj, |12| . A gaussian wave packet is defined 
by 




2vr \ / (q-Q 



,2 



4cj 2 



where a pP . is the width of the wave packet in momentum space. It is also useful to define 
the width & x p in in configuration space by a pP . a x p- = 1/2. The normalization of the states is 
defined as in Ref. ||, so that / ^iVKq)! 2 = 1- 

If the wave packet is sharply peaked around its average momentum Q, i.e. a Pm <C E P , the 
energy can be approximated by 

^„(q) = ^ in (Q) + vp in -(q-Q), 



where Ep in (q) = yq 2 + m P in an d v p 4 „ = Q/E Pin (Q). This approximation means that disper- 
sion is neglected in the external wave packets. The factors Mp(q, k) and M£>(q' , k') multiplying 
the exponential vary slowly over the width of the wave packet and can be approximated by 
their value at the average momentum. They can be factorized outside the sum over the mass 
eigenstates since we have assumed that the neutrinos are either relativistic or nearly degenerate 
in mass. 

With these approximations, the momentum integrations in the overlap function (^) can be 
done analytically. One then performs the integrations on x, x' and t, t', which are all gaussian. 
The final result can be written as 

V>(P°,P) = N exp (-/p(p°,p) - f D (p°, P )) , (6) 
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with 



Jp\P ,P) - TZo h 



4cj 2 4cj 2 ' 

^ u pp w eP 

where £p = £^ P . n — E Pout , p P = Q — K. The function f D (p°,p) is defined in the same way, 
with the index P replaced by D, except for the energy-momentum which is defined so as to be 
positive: E D = E Dout — E D . n and p D = K' — Q'. The constant N includes the normalization 
constants as well as the factors Mp y D evaluated at the average momenta. 
A new width a pP has been defined by o- pP a XP = 1/2, with 

1 1 1 

+ 



(T 2 (T 2 (T 2 

XP XP^ n "CPout 

o p p can be interpreted as the momentum uncertainty at the source. The spatial width a XP is 
mainly determined by the external particle with the smallest space width. This is expected since 
the production region depends on the overlap in space-time of the external wave packets. 
The symbol v P is defined by 

, / v P . v P . \ 

,_Z I %tl 1 r out I 



— ^ 1 in 1 

V-rP 1 



'XP \ a 2 '2 

u v> P. u 7- 1 



It can be interpreted as the velocity of the production region, approximately equal to the velocity 
of the particle with the smallest spatial spread. 

The symbol £ P , satisfying < S P < 1, is defined by 



/ v 2 v 2 

2 I Pjn 1 Pout 

P — yp I — 1- 



2 1 ^2 



XP n 



Finally, the quantity 



2 2 (\^ 2\ ^ 2 

a eP = °pp \ hp ~ V P> - a pP 



can be interpreted as the energy uncertainty at the source, or also as the inverse of the time of 
overlap of wave packets during the production process. Indeed, one can show that 

°eP = E dfer (V« " V,) 2 , (7) 

where the sum is over all wave packets involved in the production process. This sum is dominated 
by the term including the two smallest wave packets in configuration space (unless their velocities 
are nearly equal). If a x \ is the smallest width and a X 2 the second smallest, one obtains 

_ l V l ~ V 2| _ 1 (8) 

eP rpoverlap ' v ' 

where Tp Verlap is defined as the duration of the production process. Thus, a eP can be interpreted 
as the energy uncertainty at the source, since it is proportional to the inverse of the time of 
overlap of the external wave packets at the source. The quantities a XD , a pD , v D , S D , a eD , 
rpoverlap similar definitions and properties. 
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Note that stationary boundary conditions are recovered by setting v P>D = and cr eP D = 0, 
with a pP>D different from zero. Moreover, it is reasonable to impose the constraints \v PtD \T P V p rlap < 
Sp jD , where S P (resp. S D ) is the size of the macroscopic region of production (resp. detection). 
Thus we shall assume the following constraints: 



' P.D\ 



< S P , D . (9) 



These bounds are very conservative, since we shall see that stationary models such as those 



found in Refs. [16, 18] are recovered by setting \v P>D \/&ep,D = 0. In the example of the Grimus- 
Stockinger model [16], an initial stationary neutron (v>. >n = 0) decays into a stationary proton 
(v P()iitiF = 0), a 'plane-wave' electron (cr XPout:el = oo) and the intermediate neutrino. At detec- 
tion, the neutrino collides with a stationary electron (v^,. = 0) and the outgoing neutrino and 
electron are represented as plane waves (a XDouUV = a X D out ,ei = oo)- 

As the propagation distance is macroscopic, only processes satisfying global conservation of 
energy-momentum have a non-negligible probability of occurring. Since our aim is not to prove 
this well-known fact, we impose that 

Pp = Pd = Po and E P = E D = E Q . (10) 

This approximation allows to do expansions around po and Eq. An associated velocity can be 
defined by vq = pq/Eq. 



3 The oscillation formula, with or without stationary boundary 
conditions 

In this section, we shall explain a quick method to evaluate the flavor-changing probability (|5|) 
with the amplitude given by Eq. (|2|) and the overlap function by Eq. (||). Since the experimen- 
tal conditions are such that the propagating particle is on-shell, the main contribution to the 
transition amplitude (||) comes from the pole of the propagator. However, one has to be careful 
with the choice of the contour as the analytic continuation of most overlap functions diverges at 
infinity in the complex plane. The integration on the 3-momentum can be done with the help of 



the Grimus-Stockinger theorem [16|. Let ip{p) be a 3 times continuously differentiable function 
on R 3 such that ip itself and all its first and second derivatives decrease at least like 1/p 2 for 
|p| —* oo. Then, for a real number A > 0, 



, 3 V(p)e ip ' L l^oo 2vr 2 r- 
d 6 P , K ' 9 , . — > — —ij}{y/A\) 



iVAL 



A - p 2 + ie L 

where L = |L| and 1 = L/L. For A < 0, the integral decreases like L~ 2 . 

The remaining energy integral in the amplitude (^) can be done by a saddle-point approxi- 
mation [ 12 1 . However, it is quicker to perform first the time average in the probability (||), which 
yields a delta function, and makes one of the energy integrations trivial. At this stage, one has 

J dTAA* ~ -1 f dE^(E, qi 1) r(E, qj 1) e^~^ L , (11) 



8 



where ^(E,p) is the overlap function defined by Eq. @ and qj = yE 2 — m 2 . 

Eq. (|ll]) shows that the transition probability can be interpreted as an incoherent sum 
(i.e. occurring in the probability) over energy eigenstates: interference occurs only between the 
components of tp(E,p) having the same energy pq |. In this way, the correspondence between 
models with and without stationary boundary conditions is obvious: the time-integrated non- 
stationary probability is equivalent to the energy-integrated stationary probability. For example, 
the oscillation formula obtained by Grimus and Stockinger with stationary boundary conditions 
[ p!^ | has the form of the integrand in the right-hand side of Eq. (|i~l]). This equivalence confirms 
that the stationary case can be obtained from the more general nonstationary case in the limit 
of a vanishing energy width. This limit will have to be checked explicitly on the final oscillation 
formula, as it is not obvious that the approximations involved in the computations respect this 
feature. 

If the coordinate system is chosen so that L is oriented along a coordinate axis, it is easy 
to rewrite the integral (|ll]) as 

J dTAiA* ~ ^ J dE e Kn)L-f,{E)- fj {E) ^ (12) 
with the definitions fj(E) = fj P (E) + fj D (E) and 

(jE 2 -m 2 -p ) 2 (E-E -(JE 2 -m 2 - Po )v P ) 2 
fAE) = W J L + ^ W 4a2 LJ - , (13) 

where v P and po are the components of v P and po along L, while Eq has been redefined so 
as to absorb the transversal part of po • v P . The definition of fj D (E) is similar. The function 
g(l) expresses the geometrical constraint between the direction of observation 1 = L/L and the 
momentum po: 

9 (.)=exp(-(E^!2), (14) 

and restricts the neutrino propagation to a cone of axis po and angle arcsin(a p /po). The 
momentum width o~ p , defined by 

1 _ 1 1 

— g — - n I q— 5 (15) 

(T (J (T 

u p u pP u pD 

is approximately equal to the smallest width among the production and detection momentum 
widths. The associated width a x in configuration space is defined by a p a x = 1/2. 

The remaining energy integral in Eq. (12) can be performed as a gaussian integral by 
expanding the integrand to second order around the maximum of its modulus (this is called 
Laplace's method |2(J). The result takes a much simpler form if an expansion in small mass 
differences is done around (Eq, po, mo, vo), where m\ = E^—p\ and vq = pq/Eq. The expansion 
parameter is noted e and refers collectively to all 5m 2 /2EQ, where 5m 2 = m 2 — m\. We shall 
calculate the transition probability to 0(e 2 ) in the real part of the argument of the exponential 
(since the order 0(e) vanishes) and to 0(e) in the phase. The gaussian integration will be 
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consistent with the e-expansion if the extremum is computed to 0(e), the real part of the 
argument of the exponential to 0(e 2 ), the phase to 0(e), the first derivatives to 0(e) and the 
second derivatives to 0(e°). It is important to take into account all these terms, if the result is 
to coincide with the results obtained by the methods explained in the fourth section. 
The modulus of the integrand in Eq. (^) is maximal for 



E; 



Smf + 5m 2 mi „ 
E + P — 1 + 0(e 2 ' 



4E 



where the dimensionless number p is defined by 



peff 



1 



v P (v 



v D (vo 



or 



The effective width in momentum space o" pe // is defined by 



+ 



Of 



(16) 



(17) 



"peff "V 

with Op defined by Eq. (15). The effective width can be interpreted as the energy- momentum 
width of the oscillating particle, since it is the width of the overlap function. Indeed, the value 
to 0(e°) of the second derivative of fi(E) + fj(E) at the extremum reads 



2 dE 2 



\E„ 



i.i j 



(18) 



The effective width is dominated by the smallest among the energy uncertainties (recall that 



O, 



eP.D 



< a 



pP.D 



). The effective width cr X eff in configuration space, defined by o pe ffO, 



xeff 



1/2, is then approximately equal either to the production or to the detection time uncertainty, 
depending on which one is the largest. In the stationary limit (v PtD = and a eP]D = 0), cr pe ff 
goes to zero. 

The parameter p has been defined so as to be in correspondence with the notation of 



peff 



Ref. [ 12 1 . The symbol u) appearing in that article is related to our notation by u 
Note that the authors of Ref. 12] do not compute p explicitly and also take the relativistic limit 
uo = l- The explicit value of p is very interesting to know, since p = in the case of stationary 
boundary conditions, in which case all mass eigenstates have the same energy Eq. 
The value to 0(e 2 ) of fi(E) + fj(E) at the extremum reads 



2\2 



fi{Eij) + fj{E l3J 



(5mfY + (5m 2 ) 



+ 2?r 



2 I P a xeff 

T OSC 



where the mass width a m is defined by 



2 

a peff 



1 



1 



-S -5" + 



1 



Of 



or 



+ 



o 2 <J 2 



The expansion of the phase present in Eq. (12) around the extremum reads 

5m%L 



(qi -qj)L 



2p 



+ 



2p 2 v 



(E - E. 



> j i 



(19) 



(20) 



(21) 
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The second derivative of the phase is of 0(e) and can be neglected with respect to the second 
derivative of fi + fj . 

The gaussian integration in Eq. (|l^) may now be performed with the help of Eqs. (|l8|), 
( |l9|) and (21). As a result, the flavor-mixing transition probability (||) corresponding to the 



propagation over a distance L = L 1 of a neutrino, initially of flavor a, but detected with a flavor 
(3, reads 



V a ^(Ll) ~ -L exp (- (P ° 2 *> 1)2 ) E ^ ^ % 




x exp I _ 2 7ri— - I I - 2vr 2 P -^L | - l - n \ J ' I . (22) 

The oscillation length L°^ c for the masses rrij and rrij is given by 

rose _ ^Po {0 o\ 

Li i ~ tort' (23) 
where Smfj = mf — m 2 is taken to be positive. The coherence length Lf^ h is defined by 



{dm}) 2 + (5m 2 ^ 2 



r coh 1 PO rose tnA\ 

V2TT (Tpeff 



Let us now proceed to the analysis of the oscillation formula (22). The first term in the 
second exponential of Eq. ( |22| ) is the standard oscillation phase proportional to the propagation 
distance. The second term, or coherence-length term, leads to the vanishing of oscillations beyond 
the coherence length L^ h . This phenomenon is partly due to the progressive separation of mass- 



eigenstates wave packets propagating in space (see discussion after Eq. (p8|)). The third term 
is a localization term, i.e. a constraint imposing that oscillations vanish unless the oscillation 
length is larger than the space-time uncertainty: 

f'u'i />«■■< If- ( 25 ) 

This condition can be rewritten as 5m 2 -/p^ < a pe ff/p, stating that oscillations vanish if the 
energy-momentum measurements allow to distinguish between the different mass eigenstates. 
The coherence length and localization conditions were already predicted in intermediate wave 
packet models of neutrino oscillations @, fflj} . 

The fourth term in the second exponential of Eq. (p2^) could be a matter of concern since 
it does not vanish in the limit rrii = rrij. First of all, note that this kind of term is not specific 



to our computation. For example, it would appear in the oscillation formula (26) of Ref. [12] if 
the terms S a {E a ) present in Eq. (22) of that article had been expanded beyond zeroth order in 
m 2 /E 2 . It is of interest to rewrite the term under discussion as 

(5m 2 f + (8m 2 ) 2 = (Srrfyf (5m 2 + 5m 2 ) 2 

16a 2 m E 2 32a 2 m E 2 + 32a 2 m E 2 " [ > 
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The first term in Eq. ([26} ) is recognized as a localization constraint and can be rewritten as 



L^>g x , (27) 



as a m ~ v$a p whether the stationary limit is taken or not. The second term in Eq. (|2q) imposes 
that 

\mf + m? — 2rr3\ 

— 4 < *m ■ (28) 

As mo is related to the average in- and outgoing momentum by rag = Eq — Pq, condition (28) 



means that the mass eigenstates have to be on-shell with respect to (Eq,pq) within the uncer- 
tainty a m . For example, this constraint is impossible to satisfy in the mixing of relativistic and 
nonrelativistic neutrinos if the available energy-momentum is such that only the lightest neutrino 
can be produced. However it has no effect on the oscillations in the two cases considered in the 
present article, namely relativistic neutrinos or nearly mass degenerate neutrinos. Condition 
( p8|) should simply be considered as expressing the conservation of energy-momentum. Such 
kinematical constraints are usually not included in the oscillation formula, though they rightly 
belong to it. A complete computation of the transition probability should not only include this 
energy-momentum constraint, but also similar terms arising from the interaction amplitudes 
Mp^d, from the prefactor resulting from the gaussian integration and from the spin structure of 
the propagator. Actually the neutrino masses should be expected to appear not only through 
mass differences but also through their absolute values. 



The oscillation formula ( P2|) is similar to Eq. (26) of Ref. [12]. However, three interesting 
new elements appear in Eq. (p2|): the parameter p has been computed explicitly, the oscillation 
formula (^2|) is valid for relativistic or nonrelativistic (but nearly mass degenerate) neutrinos, and 
the dispersion of the 'wave packet' corresponding to the neutrino has been taken into account. 
This last fact is not at all obvious from the above computation, in which the amplitude Aj 
has been inserted into the probability before doing the integration over the energy. Had this 
energy integration been performed first, the explicit dependence of Aj on (T, L) would have been 
obtained, i.e. a gaussian spreading with T in the direction L. In section |4 . 3| this spreading will be 
interpreted as the longitudinal dispersion of the 'wave packet' associated with the propagating 
neutrino. We shall also see that computations neglecting this dispersion lead to a result differing 



from Eq. (22) in the nonrelativistic limit (see Eq. (p7|)). In contrast to the method of section 
|4.3| , the question of dispersion does not even arise in the derivation of Eqs. (pd|)-(p^) in which it 
is automatically taken into account. 

Another feature of the oscillation formula ( |22| ) is the increase of the coherence length when 
a long coherent measurement in time is performed at the detector, even if the neutrino 'wave 
packets' have separated spatially |], [l4|. In that case, the energy uncertainty at detection goes 
to zero, cr eD — > 0, so that the effective width also goes to zero, fpe// — * 0, and the coherence 
length becomes infinite, L^° h — > oo. At first sight, oscillations seem to vanish in that limit, 



because the localization term depending on p<J X eff seems to diverge |12[| . If it were true, it 
would be impossible to increase without limit the coherence length by performing long coherent 
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measurements. Note that this would be in contradiction with stationary boundary condition 
models, which have a zero energy uncertainty but an infinite coherence length. This apparent 
contradiction can be cleared up by examining carefully the term pcr xe ff in the limit a eD — * 0. 
With the help of Eqs. (|), (||) and (0), one sees that 

Hcr eD ^0 | I'd | a 
Vxeff ► < S D , 

where S D was defined as the size of the macroscopic detection region. Thus the localization 
term does not give a stronger constraint than S D < L°? c . This constraint is always satisfied, 
as it is equivalent to the constraint obtained by averaging the transition probability over the 
production region. Therefore, the coherence length can be increased without bound by more 
accurate energy measurements, contrary to what was claimed in Refs. [jl2] , |l0|| . Note that this 
is not true if the accuracy of the 3-momentum measurements is increased, as the localization 
makes the oscillations vanish when the corresponding spatial uncertainty becomes larger than the 
oscillation length. Of course the opposite conclusions would be reached if experiments measured 
time, not distance. 

As a corollary, the stationary limit of the oscillation formula ( p2|) is well-defined. Stationary 
boundary conditions, v PtD = and a eP]D = 0, lead to an infinite effective width cr xe ff and thus 
to an infinite coherence length. However, the product pcr xe ff remains finite, as we have just 



explained. For example, the Grimus-Stockinger model [16] is obtained in the limit p = 0, 
a xeff ~^ °° with pa xe ff = 0. The latter condition means that this model can be recovered 
from the external wave packet model if the stationary limit has the property v PtD /a eP ^ D = 
(see discussion after Eq. @). Though condition ( |25|) becomes ineffective in that limit, the 
localization condition ( p7|) is still present. Thus there is no contradiction between models with 
stationary boundary conditions and those with external wave packets. The former type of model 
can be obtained from the latter in a smooth limit. 

In conclusion, if the observability conditions |po X 1| < <j p , L <C Lf- h and L? sc S> a x are 
satisfied, the oscillation formula (22) reduces to the standard formula ([!]), with the additional 
property of 1/L 2 geometrical decrease. 

4 Other oscillation formulas 

4.1 The Jacob-Sachs theorem 

With result (|2~^), we have been able to reconcile the oscillation formulas proposed on the one 



hand by Giunti, Kim and Lee |11, 12 1 and on the other hand by Grimus and Stockinger [16]. More 
generally, the stationary boundary conditions models have been shown to be special cases of the 
external wave packet models. We have seen that Eq. (|22| ) reduces to the standard oscillation 
formula (jl|) if some observability conditions are verified; otherwise oscillations vanish. However, 
other authors also working in a quantum field framework have argued that formula (|22|) is not 
valid for all experimental conditions. 
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First, Ioannisian and Pilaftsis |18[] derive a formula for neutrino oscillations which exhibits 
a plane wave behaviour, if the condition L/a x <C po/<7p is satisfied (with a x <C L°!- c , as usual). 
By 'plane wave oscillations', these authors mean that the oscillation amplitude Aj does not 
decrease as the inverse of the distance, that its phase depends on the direction L as pj • L, and 
that the amplitude is not negligible in directions other than along the neutrino momentum. If it 
were true, the oscillation length measured by a specific detector would depend on the direction 
of the total momentum of the initial particles. To give a typical example, the above 'plane 
wave condition' is satisfied for atmospheric neutrinos if a x > 10~ 2 cm (with L ~ 1000 km and 
p~ lGeV). Unfortunately, this condition appears nowhere in the treatment of the previous 
section. Actually, it is not clear whether the Grimus-Stockinger theorem is valid under this 
condition since its derivation uses the large L limit. 



Secondly, Shtanov [21] argues that the relative weights of the oscillating exponentials can 
be strongly sensitive to the neutrino masses, if the source and detector are strongly localized. 
Since this author works in configuration space, the derivation of the previous section cannot be 
useful to assess his claim. 

Thus, the examination of the claims made by Ioannisian, Pilaftsis and Shtanov requires 
another method of evaluation of the amplitude (0). Instead of using the Grimus-Stockinger 
theorem to integrate on the 3-momentum, we shall first perform the integration on the energy 
with the Jacob-Sachs theorem [p4|] . This theorem is based on the assumption that the energy 
spectrum of all incident particles is limited to a finite range. Thus the overlap function ip(E, p) 
is distinct from zero only for p 2 = E 2 — p 2 within certain bounds (with E > 0). On this interval, 
ip(E,p) is taken to be infinitely differentiable. In that case, the Jacob-Sachs theorem says that 
the asymptotic value of the energy integral in Eq. (^), when T — > oo, is given by its residue at 
the pole below the real axis. Thus the evaluation of the partial amplitude @ with this theorem 
yields 



where 4>j(p) = y P 2 + m j T — p • L. For gaussian external wave packets, the overlap function is 
given by Eq. (||). In principle, this function should be cut off outside the energy range determined 
by experimental conditions so as to satisfy the conditions of the Jacob-Sachs theorem. However 
these corrections are very small and will be neglected in the computations. 

It is interesting to note that the amplitude (^) is mathematically equivalent to the ampli- 
tude obtained in the intermediate wave packet model |7j , in which the neutrino mass eigenstates 
are directly represented by wave packets. The overlap function ip(^Jp 2 + m 2 ,p) corresponds to 
the wave function of the jth mass eigenstate. Thus, it makes sense, in an external wave packet 
model, to talk about mass eigenstate 'wave packets' associated with the propagating neutrino. 
Note however that the overlap function takes into account not only the properties of the source, 
but also of the detector, which is a bit unusual for a wave packet interpretation. 
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4.2 No-dispersion regime 

4.2.1 Choice of integration method 

The integration over the 3-momentum in Eq. ( |29|) cannot be done analytically. Resorting to 
the explicit form (||) of the overlap function valid for gaussian external wave packets, we see 
that the integral ( p^ ) can be approximated by means of an asymptotic expansion for which two 
kinds of large parameters can be used. On the one hand, &pp tD and o~^ PD are large parameters 
appearing in the overlap function (|6|). They suggest a second order expansion of the integrand 
around the maximum p^- of the overlap function, followed by a gaussian integration: this is called 
Laplace's method [26|. On the other hand, T and L are large parameters appearing in the phase. 



They suggest a second order expansion of the integrand around the stationary point p c ij of the 
phase, followed by a gaussian integration: this is called the method of stationary phase pfj. The 
competition between these two asymptotic behaviors implies a detailed study of the oscillation 
of the phase around the average momentum pj . The expansion of the phase in Eq. ( ]29| ) should 
be compared with the expansion of the overlap function. Although both methods are expected 
to lead roughly to the same answer in the case of the propagation of a single particle, it should 
be checked whether the delicate compensation mechanism resulting in the oscillation phase is 
independent of the method chosen. 

The study of the overlap function amounts to the study of the argument of the exponential, 
i.e. of the function (f P + f D )(p) (see Eq. (||) with p° = ^Jp 2 + m 2 ). Using as in the previous 
section an expansion in small mass differences e = dm 2 j E 2 , the value pj minimizing f P + f D is 
given to 0(e) by 

5m 2 

Pj = Po + (au P + /3u D ) —f- , (30) 



where 



vo 

Up n 



2o e 



The associated energy Ej = y Pj + Trij can be expanded to 0(e) as 

5m 2 

E j =E + p—t, (31) 
lb® 

where p = 1 + avo • u P + (3vq ■ u D . The values of the dimensionless coefficients a and j3 can be 
computed but we shall not need their explicit expressions. It is sufficient to know that p — > in 
the stationary limit. A velocity Vj = Pj/Ej is also defined for future use. 

We are now going to approximate the overlap function as a gaussian and compute its three 
characteristic widths. At the extremum pj, the Hessian matrix of f P + f D reads to (D(e ) 



v „„ = l d 2 (f P + f D ) 

^ - 2 dp a d P b [Pj> 

5- 



2 + U p Up + u D u D , 
p 
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where a p is defined by Eq. (|15|). The matrix E determines the range of p values for which the 
overlap function ifi(yjp 2 + rn 2 , p) is not negligible. As S afc is symmetric, it can be diagonalized 
by an orthogonal coordinate transformation. The eigenvalues of S afe are 



1 

1 l/o ox . 1 



al ± = ^ + -(4 + u2)±-^( u 2 +u2) 2 -4(upX Ul3 ) 2 . 

The eigenvector associated with a 2 is in the direction of u P x u D , while the eigenvectors associated 
with a 2 ± belong to the plane defined by u P and u D . In the limit |u P | S> |u D | (resp. |u P | <C |u D |), 
the eigenvalues a 2 and a 2 _ become degenerate and the eigenvector associated with a x+ becomes 
aligned with u P (resp. u D ). This is also the case in the limit of parallel u P and u D . These limits 
are relevant to the case of stationary boundary conditions which are examined below. 

Let us choose coordinate axes (e x ,e y ,e z ) coinciding with the normalized eigenvectors asso- 
ciated with (a 2 ,a 2 _,a x+ ) respectively. The quantities (a 2 a 2 , cr 2 + ) (with a p ±a x ± = 1/2) can 
be interpreted as the momentum widths of the overlap function, since they give constraints on 
the range of p values for which the overlap function is non-negligible: 



\p x 


-Pl\ 


< 


a p 




-p)\ 


< 


CJp_ 


\p z 




< 





(32) 

The case of the stationary limit is of special interest. Recall that stationary boundary 
conditions are obtained in the external wave packet model by taking v PJ3 — > and cr ePjD — > 
with |v PD | < a ePiD S P ^ D (see Eq. @)). In this limit, the axis e 2 becomes aligned with vq, 



v x ' y ~ v%% - , (33) 
and two eigenvalues become degenerate while the third diverges: 

2 _ 2 

IE— "X ' 

4+ ^2+4 + 4- (34) 

In other words, the transversal widths (i.e. in the directions orthogonal to po) are given in the 
stationary limit by a p , while the longitudinal width (i.e. in the direction of po) is given by 
a p+ <C a p . Note that the asymptotic value of a x+ is very similar to the definition ( |i~7| ) of & xe ff 
The expansion of the overlap function has to be compared with the expansion of the phase 
4>j(p) around Pj, which reads 

Hp) = 0i(Pi) + (v,-t - L)( P - Pi ) + ^{p a - Pj)R ab (p b - p)) , (35) 

where R ab = 5 ab — VqVq. The second derivatives have been evaluated to O(e ). 
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Laplace's method will be preferred to the method of stationary phase if the phase (j>j(jp) 
varies slowly over the 'bump' of the overlap function. As a p > o p - > <r p +, the variation of the 
phase will become important in the first place along the direction p x , then in the direction p y and 
finally in the direction p z . The criterion for the use of Laplace's method in all three directions 
p x >y< z w iii thus be determined by considering the largest momentum width a p . The insertion of 
the constraints ( |32| ) into the phase (|35|) yields first order conditions for a slowly varying phase, 

\v*T- 
\v]T- 

\v*T-L z \a p+ < 1, (36) 





< 


1, 


Op- 


< 


1, 




< 


1, 



as well second order condition, 



■jjr^Sl. ( 37 ) 
Eq 



where we have used the property q a R ab q b < q 2 and the fact that a p is the largest width. For a 
given T satisfying Eq. ((3?]), it is always possible to find a range of L values so that conditions 
( |36"1) are satisfied. For other L values, the amplitude is negligible as will be checked on the result 
(see Eq. (|42|)). Thus the criterion allowing to choose the integration method is determined by 
Eq. (|37|): the integration on p will be done by Laplace's method if T < Eq/g 2 or, equivalently 
(with L = |L| and po = |po|), if 

L<^ 2 , (38) 
°p 

since conditions (|36|) impose the relation L = vqT as long as a x+ -C L. 

For <7 X + > L (stationary limit), we now show that condition (|38|) is directly obtained without 
going through condition (|37|). The overlap function imposes in that limit that |p| = \J — m 2 -, 
so that we are left with an angular integration with the angular part of the integrand given by 

exp + f P " L ) ■ ( 39 ) 

Condition (35) shows that the angular variation of the phase in p^ ) is slow with respect to the 
angular variation of the overlap function, in which case Laplace's method will give good results. 
Therefore condition (|38|) constitutes a good criterion for the use of Laplace's method whether 
the stationary limit is taken or not. 

Condition (|3g|) is usually not verified in oscillation experiments, because L/a x S> Po/cp m 
most cases. This condition is the same than the one under which Ioannisian and Pilaftsis UTi 
obtain 'plane wave' oscillations. 

4.2.2 Amplitude 



If condition (|3q) is satisfied, the evaluation of the integral (|29[) can be done by Laplace's method 
and yields 



5m 2 ^ 2 



Aj ~ a p a p .a p+ exp ( —iEjT + i Pj ■ L - I j - F j (T) | , (40) 
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where the third term in the exponential comes from (f P + /u)(pj); the parameter a m can be 
computed and has the dimension of a momentum width. The function Fj(T) is defined by 



1, _ _, t /_ . T 



F 3 {T) = -(yjT - Lf (£ + i—R j (v,T - L) , (41) 

where E ab and it! a6 are considered as matrices. In the framework of the wave packet interpretation 
developed after Eq. (29), the function exp(— Fj(T)) plays the part of the space-time envelope of 



the wave packet associated with the jth neutrino mass eigenstate. The elements of the matrix 
ReiYi + i-^^R^ 1 constrain the extent of the wave packet envelope in space-time. As T increases, 
the wave packet spreads because of the i-S-R term. Thus the dispersion of the wave packet is 



due to the second order term in the expansion of the phase 4>j(pj). Therefore, condition (37) or, 
equivalently, condition (^) means that dispersion has not yet begun in any direction, transversal 
or longitudinal. For that reason, the range of L values defined by L <Po/cri will be called the 
no- dispersion regime. Of course, this interpretation is not valid for o~ x+ > L, in which case the 
propagation time T becomes indeterminate and dispersion loses its meaning. 

Now that the origin of dispersion has been clarified, the term in Eq. ( |4~i"| ) including R can be 
neglected with respect to S. Moreover, we approximate vj by vo in Fj(T). In comparison with 
the oscillation formula (22) derived in section |3L this approximation will lead to the absence of 
the coherence-length term, since this term exclusively arises, when the dispersion is neglected, 
from the velocity difference v, — Vj. Dropping the index j, the wave packet envelope in Eq. ( ^0[ ) 
can then be written in the coordinate system diagonalizing £ as 

(vnT - L x ) 2 (ylT-Lyf «T - L z ) 2 
F ^= la* + V + \o-\ ' (42) 

^ u x ^ u x— ^ u x+ 



which shows that the conditions (pq) assumed for Laplace's method are required to obtain a 
non-negligible amplitude Aj. If cr x + "C L, condition ( |42"| ) imposes the relation 



v T = L + 0(a x+ ) . 
In that case, the phase of the amplitude (HQ) is given to 0(e) by 



Sm 2 

EjT - p j ■ L = E T - po • L + — -L {L + 0{a x+ )) , (43) 

2p 

where the expansions ( p0| ) and ( |3TD have been used. This result shows that the phase difference 
between the amplitude Ai and A* is equal to the standard oscillation phase present in Eq. ([[]). 

What is the form of the amplitude ( ^0[) in the case of stationary boundary conditions 
{°~x+ ^ L)7 This question is important since such conditions were assumed by Ioannisian and 
Pilaftsis in their derivation of 'plane wave oscillations' Jl8| . Let us first examine the wave packet 
envelope F(T). In the stationary limit (see Eq. (|34"|)), a x + becomes larger than L z so that the z 
component in Eq. (|42|) imposes that \vq\T < c x+ . Putting this condition together with Eqs. (|9|) 
and (HU), we obtain 

I x >y i 

K \ T ^ ~rw a *+ ~ Sp - D ■ 



18 



Thus the wave packet envelope (|42|) gives the following constraint in the stationary limit: 



\L x >y\<S P , D , (44) 

which means that, for propagation distances much larger than the size of the source and detector 
(L S> S P _ D ), the amplitude (^) is negligible unless the direction of observation L is nearly parallel 
to the average neutrino momentum po . Thus this constraint is valid whether the stationary limit 
is taken or not. 

Let us now examine the phase of the amplitude ( |40| ) in the stationary limit. The expansions 
(|30l) and (|3l|) and the property p — > show that 



Ej — > Eq 



5m 2 



v • Pj -> v • po - 

Using L = voL/|vo| (with correction terms given by Eq. (|44l) ), the phase of the amplitude is 
given to 0(e) by 

5m 2 - 

EjT - Pj • L = E T - po • L + — — i (L + 0(5 PiD )) , (45) 

which leads again to the standard oscillation phase present in Eq. ([[]) since L 3> S PD . Eqs. ( [44] ) 
and (|45|) also show that apart from a factor exp(-iEoT) which can be factorized from the sum 
over the mass eigenstates, the amplitude (|40[) is independent of time in the stationary limit. 

We have thus demonstrated that the amplitude (|4G|), derived under condition (p8|), is neg- 
ligible in directions other than the neutrino average momentum whether the stationary limit is 
taken or not. Therefore the oscillation phase has the standard form given in Eq. © and no 
'plane wave oscillations' can be observed, contrary to the claim made by Ioannisian and Pilaftsis 
[jig] . Also, the absence of the 1/T (or 1/L) factor in Eq. fliOD, as noted by the same authors, is 
easily understood by noting that the absence of dispersion entails that the amplitude does not 
decrease with the distance. 

4.2.3 Probability 

As in section ||, the last step towards the oscillation formula consists in computing the time 
average of the transition probability, which is a gaussian integral on T: 

( (5m 2 ) 2 + (5m 2 ) 2 (E; — E,;) 2 \ 

dT ~ ^ (-*™ - \UV " ^™ - 2F(To> j ■ (46) 

where To is the solution of F'(Tq) = 0, while the phase ^-(Tq) is given by 

^•(T ) = (Ei - Ej)To - (pi - pj) ■ L , (47) 
The function F(T) can be rewritten as 

F(T) = , (48) 
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where 



v = a x VT, 1 v 



L = a x V^L, 
with X -1 = diag(a~ 2 ,a~*,a~+). It is then easy to compute 

To = (49) 
v o 

jpfrp \ (VO >< L) 2 , . 

T"(T ) = (51) 

We now show that results fl43])-(f45|), obtained at the level of the amplitude, can be reproduced 
by a careful analysis of Eq. (|4^). This is not completely trivial as the integration variable T 
becomes indeterminate in the stationary limit. 



The insertion of expression (pG) in Eq. (46) shows that F(Tq) plays the role of a directional 



constraint. More precisely, j dTAiA*- is non-negligible if L is nearly parallel to vo or, for 
a x+ <C L, if L is nearly parallel to vo: 

1 L = P 1 L + 0{a x+ ). (52) 
l v o| 



With the substitution (p2;), the value of To given by Eq. (49) becomes To = L/\vq\ + 0(a x+ ), 
so that the phase ( |47|) becomes 

Mn) = ^r i (L + 0(a x+ )), (53) 

which is equal to the standard oscillation phase present in Eq. (|l|) and equivalent to result ((43|). 

If a x+ ^ L, v and L z tend to zero so that F(Tq) should directly be studied as a function of 
vo and L. General conclusions for arbitrary a x ± can be drawn from the study of the quadratic 
form in (L x , L y , L z ) associated with F(Tq) = 1. This analysis shows that there is an eigenvalue 
S3 = corresponding to an eigenvector along vo- The two other eigenvalues S12 are positive 
(with S2 < si) , so that the surface F(Tq) = 1 in (L x , L y , L 2 )-space is a cylinder of elliptical 
section with an axis along vo- This geometrical picture can be interpreted as imposing that the 
components of L orthogonal to vo should be smaller than yl/s2, whereas there is no constraint 
at all on the component of L along Vo- In the stationary limit, the lengthy expressions of the 
non-zero eigenvalues become 



•51 



1 

4^1 ' 



x 



S2 _> 1 a * W (54) 
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The properties of the stationary limit, namely Eqs. (||), (p3| ) and @), lead to the bound -^/l/s2 < 
S p t Q. In the stationary limit, the components L x,s/ (which are then orthogonal to vo) should 
thus be smaller than S P}D . Thus the time-averaged probability is always negligible in directions 
other than along the average neutrino momentum po- Note that this condition reproduces the 
constraint (^4|) derived in the study of the stationary limit of the amplitude (|40|). The phase 
4>ij(To) can then be evaluated as in Eq. (45), leading to the oscillation phase 

5m? ■ 

<f> ij (T ) = —^(L + O(S P , D )), (55) 
which is again equal to the standard oscillation phase present in Eq. (jl|) and equivalent to result 



Putting together the results (|46|), (|50D, (51) and (|53|) (or (|55|)), the flavor-mixing transition 
probability (||) for a propagation distance satisfying condition (^) can be written as 

7V^(L) ~ exp ^- (V 2 ° a g- L / ^) E F - V *P Via % 

x exp — 27tz 27r — - — tt^ — , (56) 

1 \ J l^lEl 1 



where a xe ff = cr x \v \/\v \ and L° sc is given by Eq. (|2j). 

This oscillation formula is very similar to Eq. (^2|). There are three main differences, two 
of which have already been mentioned. First there is no coherence-length term in Eq. (|56|), 
which is due to the neglect of the terms beyond O(e ) in the evaluation of F(To) (see Eq. (|i"2j)). 
Secondly, the geometrical decrease in 1 jl? is lacking in Eq. (^) , which is explained by the fact 
that the dispersion of the neutrino 'wave packet' is not yet significant. Finally, the directional 
constraint present in Eq. (^) restricts the neutrino propagation to a cone of angle arcsin(a p /po), 
whereas the directional constraint in Eq. (5£) confines the propagation to a cylinder of radius 
r (with a x < r < S P , D ). This different behavior is also a result of the absence of dispersion for 

Before closing this section, it is interesting to understand why a directional constraint is 
missing in Ref. for L < po/crp 1 , as this fact explains the 'plane wave oscillations' result. At 
the end of their computations, the authors of Ref. obtain an amplitude Aj whose dominant 
term depends on exp(ipj |L'|), where |L'| = VlV 2 is the complex 'norm' of a complex vector 
L' = L — 2icr 2 po (we have translated their results in our notation through the correspondence 
qj — > pj, k — ► po, 5l 2 — > 4cr 2 , L — ► I/, I — ► L). For L <C Po/&p, the quantity |L'| can be expanded 
and the argument of the exponential reads to second order 

• it o- 2 ir^o 2 , • Po-L jjj / 2 (p -L) 2 \ 

tpj\L - 2ia x p \ = 2a x p pj + ipj — L ^ . 

po 4o-2p V Po J 

The two last terms were neglected in Ref. [^] though L a x . They lead to the following 

directional constraint: 

/ (v xL) 2 \ , . 
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where the factor 4 becomes a factor 2 when the amplitude is squared. Therefore the result of 
Ref. |18| also includes for L -C po/o"p a directional constraint forbidding plane wave oscillations. 
Note that the stationary limit assumption |v PiD |/cr ePiD = leads to s\ = S2 (see Eq. (f>|)), 



so that the constraint (|57| ) becomes equal to the one present in Eq. (56). Thus the condition 
l v F,r>|/°"ep,£) = seems generic for models with stationary boundary conditions, since it was also 
applied in the case of the Grimus-Stockinger model (see discussion after Eq. (0) and at the end 
of section |3|) . 

In conclusion, the analysis of the amplitude ( f40|) (as well as of the transition probability 
(|56"|)) derived in the external wave packet model under condition (|3^) shows that 'plane wave 
oscillations' do not exist. This negative result is confirmed by a reexamination of the formula 



derived in Ref. [18]. Besides providing a rigorous method to compute the transition probability 
(56), the external wave packet model has the advantage of associating a clear physical picture to 
the different stages of computations. For example, condition ( p8| ) can be interpreted as meaning 
that the wave packet associated with the neutrino has not yet begun to spread in any direction. 
Finally, the oscillation formula ( |56"|) reduces to the standard formula (|l|) if the observability 
conditions |vo x L| < vqu x and L°? c 3> a x are satisfied. 

4.3 Transversal- and longitudinal-dispersion regimes 



Let us now assume that condition (35) is not satisfied, i.e. L > pq/ 'a 2 . In that case, Laplace's 



method cannot be used to integrate on all three components p x ^' z in the amplitude (29) since 
dispersion becomes significant. However the spreading of the neutrino 'wave packet' is not 
identical in all directions. More specifically, the onset of dispersion in the direction po can be 
delayed by two factors. First, the matrix element R ah present in Eq. (]35| ) leads to a relativistic 
contraction (of 1 — Vq) in the direction po of the dispersion of the neutrino 'wave packet' (see 
Eq. (^j])). Secondly, the momentum width along po is given for a eP , D <C a x (i.e. in the stationary 
limit (|34|)) by a vanishing <t p +. Thus Laplace's method will be valid for a longer time T in the 
direction po than in directions transverse to this vector. For this reason, the choice of the 
integration method in the direction po will be postponed for a short while, while the method of 
stationary phase will be preferred for momentum integrations in directions transverse to po- 



Let the z axis be along L, i.e. L = Le z . As in section 4.2, the examination of the amplitude 
( p9|) shows that the quick variation of the phase averages the amplitude to zero unless po and 
L are nearly parallel. The method of stationary phase will thus be applied in directions p x ' y , 
the stationary points of which are given by p x = p y = 0. The result of the method of stationary 
phase for the transverse momenta in (^) can be written as follows: 

Aj ~ =r^- f dpe-^M-JiGO, ( 58 ) 
T -in J 

where p = p z and = Eo/2(Tp. The geometrical constraint g(l) is given by Eq. (|i~4l), while the 
momentum width a p is defined by Eq. Qi~5|). The phase (j>j(p) is given by 



4>j (P) = \/P 2 + m 2 j T-pL. (59) 
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The original overlap function is partly included in g(l) and partly in 



flip) = fjpip) + fjn(p) , (60) 

with 



p / \ (P-Poj 
fjpip) = + 



p 2 + m 2 - E - (p - po)v P 



where v P = vp and po = p$. The energy Eq has been redefined so as to absorb a factor 
PgVp + p\vp. The definition of fj D (p) is similar. As T > E^/a 2 e jp the term 1/(T — ifj,) 
can be approximated by 1/T so as to provide a prefactor 1/L 2 in the transition probability. 
This expected geometrical decrease is seen to originate in the transverse dispersion of the wave 
packet corresponding to the oscillating particle. For future use, we define a reference mass mo 
by m-0 = Eq — Pq and a velocity vo by v o = Po/Eq . 

As regards the longitudinal momentum integral, it is tempting to proceed as in section ||, i.e. 
to compute first the time average of the transition probability before integrating on p. However, 
the prefactor 1/(T — i/x) present in the amplitude (|58|) gives in that case a delta function look- 
alike of width pT x ~ (Tp/i^o, introducing an additional momentum uncertainty which is larger 
than the mass difference 5m 2 /E since L°? c > po/a 2 . For this reason, it is preferable to avoid this 
shortcut (note that it yields the same final answer as given by the following method). Moreover, 
it is interesting for the physical interpretation to postpone the time average, so as to obtain 
the explicit dependence of the amplitude on time and distance. As in the case of transverse 
momenta integrals, the choice of the method to perform the longitudinal momentum integration 
([58]) is done by comparing the expansions of the phase and of the overlap function around the 
value pj for which fj(p) is extremal. One obtains to 0(e) 

5m 2 : 

Pj = Po + (P ~ 1) , 
*Po 

where the dimensionless number p is defined by Eq. ([!(]). The associated energy Ej = ^Jp 2 + m 2 
and velocity Vj = Pj/Ej are given to 0(e) by 

5m 2 

Ej = Eo + P 2E ' 

v, = -o + (Kl-,o 2 )-l)^4- (61) 

The expansions of fj(p) and <j)jip) are given by 

f.( p ) c (P-Pi) 2 
IAP) ~ Aa 2 

m 2 T 

MP) = MPi) + ( v 3 T - L )(p-Pj) + ^3(P-Pj) 2 , ( 62 ) 
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where a pe f f is given by Eq. (|17|) and can again be interpreted as the longitudinal width of the 
overlap function. Laplace's method will be used if the phase (p^) varies slowly over the width 
(Tpeff, i-e. if the two following conditions are satisfied 

\vjT - L\2a peff < 1, (63) 
-^Heff S !■ (64) 



As in section the first order constraint ( p3|) is included in the result of Laplace's method. 
Thus the criterion allowing to choose between Laplace and stationary phase methods is given by 
Eq. (|64[). In other words, it is better to use Laplace's method if T is smaller than a dispersion 
time T- isp defined by 

rpdlSp _ ^0 

The term 'dispersion time' is justified by the fact that it is the time at which the longitudinal 
dispersion of the amplitude becomes important, more precisely twice the initial size. A dispersion 
length L^ tsp can be defined by L^ lsp = voTj isp . The distance range Po/&p e ff ^ L < L^ tsp will 
be called the tranverse- dispersion regime. For L > L^ isp , the stationary phase method is more 
accurate: this distance range will be called the longitudinal- dispersion regime. Various estimates 
of the dispersion length are given in ||, showing that the concept of dispersion length is relevant 
to nonrelativistic particles as well as to supernova neutrinos, and possibly to solar neutrinos. 

In the transverse-dispersion regime, the evaluation of the amplitude ( |58|) as a gaussian 
integral around pj gives 

g{\)o- peff (-ET+- L- ( J?iX 1 ^ T ~ L f I ((T) 

3 ~ T^l + iT/Tf sp ^ \ 1 3 W3 \^ E o ) 1 + iT/Tf sp ^l eff ' ' 

The width a m is defined by Eq. (|2^). The amplitude (|6^) behaves as a wave packet of group 
velocity vj and space-time extent (1 + (T /T^ lsp ) 2 ) 2 cr X eff- If the longitudinal dispersion is 
neglected (Tj tsp = 00), the amplitude ( plf ) is similar to Eq. (18) of Ref. (T^j. 

The following step consists in computing the time averaged transition probability (||). If 
dispersion is taken into account (rj lsp < 00), a tedious gaussian integration by Laplace's method 
(see H) yields the oscillation formula (^). If the dispersion is neg lected (Tf sp = 00), a much 
shorter computation leads to the same oscillation formula (]22|), except that the following sub- 
stitution has to be made: 

L t^ 1 n % il {FALSE) . (67) 
IPC 1 — ^0) — !| 

The incorrect multiplying factor has its origin in Eq. (pl|): in the limit Tj lsp = 00, the coherence 
length term arises only from the difference between the group velocities Vi and vj. However the 
factor \p(l — Vq) — 1| tends to 1 in the relativistic limit, so that the substitution (|67| ) becomes 
trivial. This observation explains why our result ([22] ) coincides with Eq. (26) of Ref. [p^ l, as 
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the authors of this article, while neglecting dispersion, consider only relativistic neutrinos. Note 
however that even relativistic neutrinos spread at large distances so that a calculation neglecting 

- disp 



dispersion such as in Ref. |C| is only valid for L < L- 

At sufficiently large distance, dispersion becomes significant and all neutrinos propagating 
freely enter into the longitudinal-dispersion regime. In this regime, we have argued that the 
integral (|5Sj) should be evaluated with the method of stationary phase. The stationary point of 
the phase 4>j{p) is given by 

Vcl 

Pd,j = rn 



where v c \ = L/T. It can be interpreted as the classical momentum of a particle of mass rrij, 
travelling at the classical velocity v c [. The evaluation of the amplitude ( |58| ) as a gaussian integral 
around pdj gives 



a g(l)a pe // 
Aj ~ — — exp 



TJl + iT/T* 



disp 



-" n i \ 11 11 fiM + Vpeff - ~ . T/T dis P 



i \ 1 > ) 



(68) 



where fj(p) is defined by Eq. (pQ), /j(p) refers to its derivative. The wave packet interpretation 
of the amplitude (^) is not obvious but the shape of the associated wave packet can be studied 
by an expansion around the maximum of the amplitude. Once more, the time average of the 
transition probability (|5|) is computed by Laplace's method (see Ref. ||), yielding again the 
oscillation formula (|22|). 

It is striking that the two different methods of approximation used in this section, a priori 
valid in different regimes, give the same oscillation formula, which also coincides with the result 
obtained with the approximation method used in section ^, Eq. (|22|). The dispersion length 
does not play any special role in the final result. Thus, each method is accurate enough to be 
extended to the whole range of distances. However, the physical interpretation of the terms 
appearing in the oscillation formula (p2|) depends on the relative values of L and L^ p . From 



the detailed analysis done in f6j, these terms can be interpreted as follows: 

• In the transverse-dispersion regime, the coherence-length term mainly arises from the 
progressive spatial separation of 'wave packets', due to their different group velocities. 
The localization term mainly arises from the initial spread of the 'wave packet', which 
must be smaller than the oscillation length. This is in agreement with the traditional 
explanation of these two terms. Note that ]J^ sp <^ j^ s v ^ m . ^ m ^ j n cas6) one 
can check that decoherence occurs for L < ^ lsp . 

• In the longitudinal-dispersion regime, the coherence-length term mainly arises from the 
spreading of the neutrino 'wave packet': the interference term is averaged to zero by 
the time integral as the spread of the wave packet becomes large. The localization term 
mainly arises from the separation of the 'wave packets': if they are spatially separate at 
the onset of the longitudinal-dispersion regime, they will never again overlap because their 
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separation increases as fast as their dispersion. This is a new explanation for the origin of 
these two terms. 



Thus the origins of the coherence-length and localization terms in the transverse-dispersion 
regime become respectively the origins of the localization and coherence-length terms in the 
longitudinal-dispersion regime. The situation can be summarized in the following diagram, 
where 'wp' is an abbreviation for 'wave packet': 





transverse-dispersion regime 


longitudinal-dispersion regime 


Coherence length 
Localization 


decreasing overlap of wp 
initial spread of wp 


increasing dispersion of wp 
constant overlap of wp 



TABLE I. Origins of the coherence-length and localization terms. 



It is now possible to examine Shtanov's analysis in configuration space [21|, since the explicit 
L, T dependence of the amplitude Aj has been computed for the whole range of distances 
(see Eqs. (f5]), (|66|) and (|68|)). This author works directly with the neutrino propagator in 
configuration space, and considers the space-time variables x, x' (see FIG. 1) as macroscopic 
variables. His model will be considered in a scalar version in order to make the comparison with 
our results easier. In configuration space, the propagator for a scalar neutrino propagating from 
x to x' is asymptotically (rrij^/ (x' — x) 2 3> 1) given by 

J[. ^ p-imjy/{x'-x) 2 



(69) 

_ 3)2)3/4 

where {x 1 — x) 2 is the Lorentz interval. Shtanov then computes the convolution of this propagator 
with a source. In the case of a monochromatic source, he obtains the standard oscillation 
formula. However, a convolution with a strongly localized source (& Xt t < ^/E) leads to an 
amplitude that keeps its dependence on the mass prefactor ^yrrij. In that case, the transition 



probability is not equivalent to the oscillation formula (|2^), at least if the neutrino masses are 
not nearly degenerate. We proceed to show that Shtanov's result is incorrect. Note first that 
the amplitude (|69| ) is in correspondence with our amplitude (|68|), computed with the stationary 
phase method, since the prefactor in Eq. (|68l ) can be rewritten, for T ^ T c - 

Vpeff V™j 



disp 
3 ' 



as 



1 

T 



1 + iT/Tf sp 



(T 2 - L 2 ) 3 / 4 ' 



where the definition ( p5|) of the dispersion time has been used. This prefactor coincides with the 
prefactor in Eq. (|69|) . However, the subsequent time average of the transition probability com- 
pletely cancels this dependence on the mass, yielding Eq. (^). This can be seen by expanding 
the argument of the exponential in Eq. ( |6"8| ) around the average propagation time. The width 

(where j c i is the Lorentz 



with respect to T of the amplitude Aj is found to be equal to 



factor associated with the velocity v c i 
cancels the rrij prefactor in the gaussian integration 



L/T), thereby providing a m^-dependent factor that 
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Shtanov does not obtain such a result, since he does not average the probability on time. 
Besides, the mass prefactors would remain even if a time average were performed on the prob- 
ability: no 'wave packet envelope' appears indeed in Shtanov's amplitude, as this author does 
not try to compute the explicit convolution of the propagator (^) with a wave packet source of 
arbitrary width. For that reason, he wrongly concludes that the mass prefactor remains if the 
source is strongly localized. Shtanov also derives another oscillation formula for mass eigenstates 
such that m2\J [x 1 — x) 2 C 1 < m\y (x' — x) 2 . However, it can be checked that decoherence 
occurs in that case. In conclusion, Shtanov's computations in configuration space do not lead 
to new oscillation formulas. 

The last case to be considered is the oscillation formula obtained by Blasone, Vitiello and 
collaborators (hereafter abbreviated as BV). These authors have attempted to define a Fock 
space of weak eigenstates and to derive a nonperturbative oscillation formula [22]. They define 
flavor creation and annihilation operators, satisfying canonical (anti)commutation relations, by 
means of Bogoliubov transformations. Apart from the speculative nature of the undertaking, 
the drawbacks of the approach are the dependence on time, not on space, of the oscillation 
formula, as well as the lack of observability conditions. The latter problem is important since 
it determines whether the new features of the BV oscillation formula are observable in practice. 
Since the BV oscillation formula tends to the oscillation formula (^2|) in the relativistic limit or 
if the mass eigenstates are nearly degenerate (with the coherence-length and localization terms 
removed), we can focus on the case of a nonrelativistic oscillating particle having very distinct 
mass eigenstates. In that case, p ~ 5m 2 /2E, so that either a p < 5m 2 /2E or p < a p . Under 
these conditions, the oscillation formula (EQ) shows that oscillations vanish, either because of the 
localization term or because of the coherence-length term (though Eq. (|22| ) is strictly speaking 
not valid in the case considered, localization and coherence-length conditions are generic features 
of quantum-mechanical models of oscillations |7j ) . Once the oscillation terms have been averaged 
to zero both in the BV formula and in Eq. (p^), these two oscillation formulas do not differ 
anymore. Therefore, the BV formalism does not seem to be relevant to the phenomenology of 
oscillations on macroscopic distances. This observation does not detract from the theoretical 
worth of that approach. 

We shall not treat here the charged lepton oscillations obtained by Srivastava, Widom and 
Sassaroli [E(J. The reader is referred to || for a discussion of this subject. 



5 Conclusion 

In this article, we have shown that all correct results obtained by applying quantum field theory 
to neutrino oscillations far from the source are incorporated in the oscillation formula (p2|) . 
Closer to the source (i.e. for L < pq/ct 2 ), this formula is modified in only one respect without 
repercussions on the oscillation length: the global spatial distribution of the detection probability 
has the form given in Eq. (^), with the consequence that the neutrino propagation is restricted 
within a cylinder instead of a cone. The quantum field oscillation formulas ( [22] ) and ( |56P are 
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equivalent to the standard result (||) if some observability conditions are verified. The new 
oscillations formulas proposed by Ioannisian and Pilaftsis [18| and by Shtanov have been 
disproved, while the Blasone-Vitiello formula [p2| has been shown to be phenomenologically 
equivalent to the standard result. Besides, arguments against oscillations of charged leptons, as 
proposed by Srivastava, Widom and Sassaroli [20|, can be found in Ref. ||. 

We have insisted on the fact that the flavor-mixing transition amplitude can be inter- 
preted in terms of wave packets, so that oscillations can be seen, like in the intermediate wave 
packet model, as the result of an interference between propagating wave packets. Moreover, the 
equivalence (|ll|) shows that this physical picture still holds in the case of stationary boundary 
conditions, provided that an incoherent sum over the energy is performed, contrary to what was 
claimed in Ref. [17]. As a byproduct of our analysis, the propagation of a particle has been 
shown to go through three regimes, separated by two thresholds which are marked by the onset 
of first the transversal dispersion and then the longitudinal dispersion of the associated wave 
packet. This picture has been found to be useful for the understanding of the terms present in 
the oscillation formula as well as for the computation of the correct spatial distribution of the 
probability. We believe that the wave packet interpretation is the only consistent picture of the 
oscillation process. In particular, it shows the irrelevance of discussions bearing on the neutrino 
energy-momentum, since it has not a unique value for a given mass eigenstate, but rather a 



spread described by the overlap function V"(y P 2 + m2 >P)- 

One could wonder what advantage the external wave packet model has over the intermediate 
wave packet model, apart from being the only consistent way to derive the oscillation formula. 
Actually, the quantum field formula ( |22| ) has two important new features with respect to the 
quantum mechanical result. First, the energy uncertainty is not put in by hand, but is defined 
in terms of the 3-momentum widths and velocities of the in- and outgoing wave packets (see 
Eq. (^)), with the result that the stationary limit is well-defined and that it is possible to analyze 
the dependence of the coherence length on the accuracy of energy measurements. Secondly, the 
oscillation formula is valid as well for neutrinos as for K and B mesons, putting it on a much 
firmer phenomenological basis. This might be useful in the light of the numerous nonstandard 
oscillation formulas existing in the literature (see Ref. || for a review). 

The stationary limit has been considered in our computations to show the great generality of 
the external wave packet model. The stationary assumption is not problematic for the derivation 
of the oscillation formula. However, this limiting case is rather unphysical. The computations 
of section ^ have indeed shown that boundary conditions can be considered as stationary when 
^xeff ^ L. This constraint, which can be rewritten either as a e p,D ^ l/£ or as Tp V p rlap > T 
(verified either at the source or at the detector), is extremely stringent. In the example of 
atmospheric neutrinos, the process can be considered as stationary if a eP>D < 1CP 19 MeV or 
equivalently if T° v £ rla P > io~ 3 s. These conditions are far from being realized at the microscopic 
level in any known neutrino source or detector. Although most neutrino sources are stationary 
from a macroscopic point of view, there is no reason to think that individual particles in the 
source and detector remain unperturbed in coherent states over macroscopic time scales: the 
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Sun is certainly not stationary at the atomic scale, and neither is a detection process where 
charged leptons are observed with a finite energy and time spread |Q . 

A subject deserving further investigation in the quantum field framework is the influence 
of the decay width of the source of the neutrino. The case of decay at rest has been studied in 
detail with the Wigner-Weisskopf approximation in Ref. [jl^J, while the treatment of decay in 
flight [p8l is still unsatisfactory. Unstable neutrinos have been studied in Ref. 

We would be sorry to leave the reader with the impression that the quantum field approach 
to oscillations is complicated. Actually, this is a consequence of the generality of the external 
wave packet model. For pedagogical purposes, it is interesting to resort to a simplified quantum 
field model, where the source and detector are perfectly located in space and in- and outgoing 
states are stationary . 
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